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Abstract 

Using the neutrino oscillations and neutrinoless double beta decay ex- 
perimental data we reconstructed an upper limit for the three generation 
neutrino mass matrix. We compared this matrix with the predictions of the 
minimal supersymmetric(SUSY) model with R-parity violation(^p ) and ex- 
tracted stringent limits on trilinear coupling constants Ajss , A -33 . Intro- 
ducing an additional U{l)x flavor symmetry which had been successful in 
explaining the mass hierarchy of quarks and charged leptons we were able to 
relate various ^pparameters. In this model we found a unique scenario for the 
neutrino masses and the I^p couplings compatible with the neutrino oscilla- 
tion data. Then we derived predictions for certain experimentally interesting 
observables. 



* Oliver . haug@uni-t uebingen . de 



1 



I. INTRODUCTION 



It is a common belief that the existence of neutrino oscillations point to physics beyond 
the standard model (SM). Recent Super-Kamiokande results strongly support the existence 
of neutrino oscillations by the observation of the zenith-angle dependence of the high en- 
ergy atmospheric z/^ events. Other hints for this phenomenon come from the solar neutrinos 
and the accelerator LSND [Q- ||^ neutrino oscillation experiment. 

The neutrino data were extensively used for testing various models of physics beyond 
the standard model |J7[. Recently there was a growing interest in the description of neutrino 
properties in supersymmetric models with R-parity violation (R-parity violating Minimal Su- 
persymmetric Standard Model ^pMSSM). It was realized a few years ago that the ^pMSSM 
framework is rather adequate for this purposes. A non-trivial Majorana neutrino mass ma- 
trix is a generic feature of the ^pMSSM as a consequence of the lepton violating ^pcouplings 

m- 

In the present paper we are studying the impact of neutrino oscillation data on the three 
family neutrino mass matrix and on the flavor structure of the MSSM couplings. 

There exists a controversy, whether a three neutrino family scenario is able to accom- 
modate all these data or an additional fourth light sterile neutrino must be included in 
the theory []TU[. Recently it was argued that three neutrinos are enough for a reasonable 



description |]rT|- of all the above cited neutrino oscillation data. A especially good fit to 
the data was obtained by taking out the LSND points from the analysis. This is motivated 
by the opinion that the LSND result needs an independent confirmation. 

Here we accept the three family neutrino scenario. In section |I| we start with consider- 



ation of the constraints imposed on this scenario |[Tl|- [0] by the neutrino oscillation data 
and show that in this framework one needs additional information on the overall neutrino 
mass scale in order to determine the neutrino mass matrix. Towards this end we use the 
neutrinoless double beta decay (OuPP) experimental constraints on the average neutrino 
mass (nil,) and reconstruct the three family neutrino mass matrix with the maximal entries 
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allowed by these data. 

In section |T| we consider the three family neutrino mass matrix in the I/ip MSSM [|15 



We allow a most general case of explicit -R-parity violation in the superpotential and the 



soft SUSY breaking sector |T6| , p!7| taking into account both the trilinear and the bilinear 
terms. In this model the neutrinos acquire masses at the electroweak scale via tree 
level neutrino-neutralino mixing as well as via 1-loop corrections We compare the 

total 1-loop three neutrino family mass matrix of the Ijip MSSM with the maximal mass 
matrix derived in section |T| and extract stringent constraints on the trilinear ^p couplings 
Aj33, A^33. Similar constraints on the bilinear parameters and products of certain trilinear 



-^p couplings from the neutrino oscillation data were previously derived in ref. |jT8[- pO |. 

It is known that the predictive power of the ^pMSSM is quite weak due to the presence of 
many free parameters. In section |^ we consider a model based on the presently popular idea 
of the horizontal U{l)x flavor symmetry. Being imposed on the ^p MSSM this symmetry 
relates many parameters and allows one to very successfully describe the quark and the 
charged lepton masses and mixing angles pT| , p2| . The Ijip couplings are also subject to 
U{l)x symmetry relations. These relations restrict the ^pMSSM neutrino mass matrix so 
that the overall neutrino mass scale becomes fixed only by the neutrino oscillation data. 
We find a unique solution for the neutrino masses and the trilinear Ijip couplings for every 
oscillation analysis. On this basis we predict the average neutrino mass {my) and shortly 
discuss prospects for the future Oz//?/3-decay experiments. 



II. NEUTRINO OSCILLATIONS AND NEUTRINO MASS MATRIX. 
PHENOMENOLOGICAL TREATMENT 

Neutrino oscillations can occur if neutrinos have a non vanishing rest mass and their weak 
eigenstates |z/°), a = e,fi,T, do not coincide with the mass eigenstates |z/j), i = 1,2,3. The 
unitary mixing matrix U that relates the weak and mass eigenstates can be parameterized 
in the three family scenario by the three angles 612, ^13, ^23- Assuming that CP- violation is 
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negligible one gets 



C12C13 



■5l2Cl3 



u 



-S12C23 - C12S23S13 C12C23 - S12S23S13 S23C13 



(1) 




where Sij and stand for sin6'jj and cos 6ij respectively. Then a neutrino produced in the 
weak eigenstate |z/°) changes its flavor content when propagating in space. The probability 
of finding a neutrino produced in the flavor state a at a given distance L from the production 
point with the energy E in the flavor state f3 is given by 



Here Amf^- = |mf — m|| is the difference of the squared masses of the neutrino mass eigen- 
states i and j. The phenomenological neutrino mass matrix Ai^'^ in the flavor space is 
connected to the physical neutrino masses rrii by the mixing matrix U as follows 



Using eq. one can extract from the oscillation experiments the mixing angles 6ij and 
the squared mass differences Amjy This information is not sufficient for the restoration 
of the neutrino mass matrix Ai^^. The overall mass scale as well as the CP eigenvalues 
(q^p of the neutrino mass eigenstates (+1 or -1 in our model) remain undetermined. To 
fix this ambiguities one needs additional experimental information other than the neutrino 
oscillation data. This information exists in a form of upper limits on the neutrino masses or 
their combinations. Experiments, measuring the neutrino masses directly, offer at present 
too weak limits, leaving the neutrino mass matrix very uncertain. Much better result can 
be achieved using the neutrinoless double beta decay (OuPP) constraints on the generation 
average Majorana electron neutrino mass 




i<j=i 





(3) 




(4) 
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From the currently best experimental limit on OuPP-decaj half-life of Ge |23 



T°/f'^(0+^0+) > 1.1 X 10^5 years (90% C.L.) one obtains (m^) < 0.62 eV [^,0. 

Now with this additional input limit we can find the maximal allowed values for the 
matrix elements m^"^ of the neutrino mass matrix. In our numerical analysis we are search- 
ing for these maximal values over the whole allowed mass parameter space. In doing this 
we take care of all the possible CP-phases of the neutrino mass eigenstates. The resulting 
absolute values of the matrix elements of this "maximal" neutrino mass matrix are 

^ .60 .97 .85 ^ 

\m'--\ = .97 .76 .80 eV. (5) 
.85 .80 1.17^ 

Here we used the results of the phenomenological analysis of the neutrino oscillation data 
(including the LSND data) made in refs. [|lT|- [|l^]. In eq. (^ the worst case of the weakest 
bounds is given. This "maximal" neutrino mass matrix can be used to test various theoretical 
approaches and allows one to extract limits on certain fundamental parameters. Below we 
are studying in this respect the ^pMSSM and find new limits on the parameters. 
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III. NEUTRINO MASSES IN MSSM 

The MSSM is the minimal supersymmetric extension of the SM. Assuming that R-parity, 
defined as Rp = (— i)3-B+l+2S ^j^^ g^j^^^ g ^j^g baryon, lepton numbers and the spin), is 
conserved one ends up with the superpotential 

Wr^ = Xf,H,L,E^ + Xf^H.Q.D'; + X^H^QiU'r + ^^H^H2. (6) 

Here L, Q stand for lepton and quark doublet left-handed superfields while E'^, f/^, for 
lepton and up, down quark singlet superfields; Hi and H2 are the Higgs doublet superfields 
with a weak hypercharge F = — 1, +1, respectively. 

In the MSSM with conserved R-parity neutrinos remain massless particles as in the 
SM. This follows from the fact that in this framework there is no room for gauge invariant 
neutrino mass terms of either Dirac or Majorana type. 
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Since R-parity conservation has no robust theoretical motivation one may accept an 
extended framework of the MSSM with R-parity non-conservation (^pMSSM). In this case 
the superpotential W acquires additional terms ||T5| 



W^^ = \^,kUL,El + K.kUQjDl + ^^,L,H, + K.i^'rDpi (7) 

so that the MSSM superpotential isW = Wr^ + Wjfi^. 

In general, R-parity is also broken in the "soft" SUSY breaking sector by the scalar 
potential terms 

F^f = K,uUL,El + KikUQ.Dl + K,,Utbpl + fxirL.H, + -^iL.Hl + H.c. (8) 

The terms in (|^ and @ break lepton and baryon number conservation. The tilde indicates 
that one includes only supersymmetric fields. To prevent fast proton decay one may assume 
A" = A" = that is commonly expected as a consequence of certain symmetry like baryon 
parity 



The R-parity conserving part of the soft SUSY breaking sector includes the scalar field 
interactions 



V^°J'= ml\(t)i\^ + \''A''LH^LE'' + X''A''HiQb''+ (9) 

=scalars 

+X^A^H2QU'' + HBH1H2 + H.c. 



and the "soft" gaugino mass terms 

1 



^GM — "2 



M^BB + M^W'^W'' + M^g^g" - H.c. (10) 



Here Ms^2,i stand for the "soft" masses of the SU{3) x SU(2) x f/(l) gauginos g, W, B while 
rrii denote the masses of the scalar fields. 

In the above sketched framework of the MSSM neutrinos are, in general, massive. 
Generically, one can distinguish the following three contributions to the neutrino masses: 

• Tree-level contribution: 

The bilinear terms in eqs. (0), lead to terms in the scalar potential linear in 
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the sneutrino fields and thereby a non-vanishing vacuum expectation values(VEVs) of 
these fields (z>j) 7^ 0. This leads to a mass term for the neutrinos by mixing with the 
gaugino fields and W^. The term fiiLiH2 in equation gives an additional mass 
term from the mixing of neutrinos with the neutral Higgsino fields H^, H2. The so 
generated non-trivial 7x7 mass matrix in the basis {u^, z/^, i^r, B^, H^, H2) can be 



brought into a block diagonal form and an effective neutrino mass matrix at tree 
level M**"*^^ can be extracted. In leading order in implicitly small expansion parameters 
{i>i) /M SUSY, IJ'i/M SUSY one gets the expression p7| : 



M*7 = ZiA«A^, A„ = /i(z>„) - (11) 

Ml + tan^ ewM2 



4(sin2/5 ^ (Mi + tan^ dwM2) - Mi M2 /x^) 



Here c/2 is the SU{2) gauge coupling and tan/5 = {Hl)/{H^). The "soft" SUSY 
breaking gaugino masses M2,i and the superpotential parameter fi are usually assumed 
to be not too far from the characteristic SUSY breaking scale Msusy ~ lOOGeV. 

gg-loop contribution: 

Another contribution to the neutrino masses arises due to the quark-squark self-energy 
loops, coming from the term \[-^LiQjD'f^ in the R-parity violating superpotential in 
eq. (0). The corresponding diagram is shown in fig. 1(a) and its contribution to the 
neutrino mass matrix is given by 



_ V- ^Kki>^'jmn MtMii{A^i + /xtan/j) 

F 

= ^2A-33A'-33. (12) 



— 6ti^ m. 



3A:33A;33mg(A^ + /itan/j) _ ^ ^, 

= ^2A,33A,33- 



In the numerical analysis we assumed the down quark mass matrix M'^ to be diagonal 
and keep only the dominant contribution of the 66-loop which is proportional to m^. 
The dominance of the heaviest internal fermion line holds for non-hierarchical X'^^f^ and 
Afj in a sense that they do not strongly grow with increasing generation indices. We 



also assume these quantities to be real. The factor 3 appears from summation over 
the (s)quark colors in the loop. 

//-loop contribution: 

The last contribution to the neutrino mass matrix at the 1-loop level is given by the 
diagram in fig. 1(b). It is induced by the LiLjE^ term in eq. (|^. It has the same 
structure as the above discussed gg-loop. The //-loop contribution to the neutrino 
mass matrix reads 

j^il ^ KkiXjmn Ml^M-^iiAli + /itan/?) 



^ 87r2 ml 

k,l,m 

A^33AJ33"^^(^f + /itan/3) _ 

^ ^2 = ^3Ai33Aj33- (i-^j 

Here, as in the case of the gg-loop, we assumed absence of the generation index hier- 
archy in Ajjfc and Afj and the charged lepton mass matrix M'^ to be diagonal. For the 
numerical analysis we kept only the dominant rf-loop contribution. 

The quantities and in eqs. (|12D, (0) denote the left-right averaged square of 
squark and slepton masses respectively. 

Summarizing, we write down the 1-loop level neutrino mass matrix TW*^ in the ^pMSSM 

as 

From eqs. ([TT|)-(p!3D one sees that this matrix is of the from 

3 

^• = E«'4' (15) 

k=l 

with the three terms in the sum corresponding to the three terms in eq. ( |I4D built of the 
three different 3-dimensional vectors a^. A matrix with such a structure has no eigenvector 
x^'^^ with eigenvalue zero and therefore all the three neutrinos have non-zero masses. This 
follows from the fact that the zero-eigenmass condition ■ f = requires the vector 
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x*^") to be simultaneously orthogonal to the three different vectors a!' {k = 1, 2, 3) which is 
impossible in the 3-dimensional space. The same arguments show that neglecting one or two 
of the terms in eqs. (|14|), (|15D results in one or two zero-mass neutrino states respectively. 



Thus in order to keep all the neutrinos massive we retain all three terms in eq. (0). 

In section || we extracted the limits on the matrix elements of the neutrino mass matrix 
in the three family scenario. Now we can translate these limits to limits on the trilinear 
couplings of the -^pMSSM. Towards this end it is enough to analyze only the diagonal matrix 



elements. Note that in case when all the three normalization factors ^1,2,3 in eqs. ([Tl])-(|T3|) 
have the same sign no compensations occur between different terms contributing to these 
matrix elements since they would be a sum of positive(negative) terms. This would allow 
one to get an upper bound not only for the sum but also for each term of the sum separately. 
As follows from the Renormalization Group Equation analysis this condition is satisfied 
for a quite wide region of the MSSM parameter space but not everywhere. In our order of 
magnitude analysis it is enough to assume that there are no large compensations between 
the different terms. For the MSSM parameters we take A ~ /x ~ rrii ^ ruf ~ Msusy and 
tan (3 = 1. The characteristic SUSY breaking mass scale Msusy is usually assumed to vary in 
the interval 100 GeV < Msusy < 1 TeV motivated by non-observation of the superparticles 
and by the "naturalness" arguments. With this choice of parameters we obtained the upper 
bounds for the trilinear couplings shown in table |. There we also display the existing 
bounds for this parameters . One sees from the table | that our limits are more stringent 



than the previously known ones. 

IV. /^pMSSM WITH FAMILY DEPENDENT U{l)x SYMMETRY 

Now suppose that the I/ip MSSM Lagrangian is invariant under the family dependent 
U{l)x symmetry. Presently this is a popular idea which allows one to predict the hierarchical 
structure of the charged fermion mass matrices and the fermion masses are in agreement with 



low-energy phenomenology pl|j2^ . In this ^pMSSMx[/(l)x-model the trilinear ^pcouplings 
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are forbidden by the U{l)x symmetry and appear when it is spontaneously broken. These 
couphngs can be generated by the effective operators 

existing in the U{l)x symmetric phase and originating from the underlying theory at the 
large scale Mx- In eq. (|16D we denoted the U{l)x charges of the quark and lepton fields 
as li,qj,dk,ek. The SM singlet field 9 has the U{l)x charge -1 and, acquiring the vacuum 
expectation value (6) ^ 0, breaks the f/(l)x symmetry. As a result the effective operators 
in eq. (|T6|) generate the following effective couplings p!7|j30|] 



A,,fc ~ e'-'°Af„ (17) 
- e'-'nf,. (18) 

Here is e = {6)/Mq and U = \li + with li,h2 (i=l,2,3) being the U{l)x charges of the 
lepton and H2 Higgs fields. The parameter e ~ 0.23 and the relative U{l)x charges 

I/1-/3I =4, 1/2 -/"si = 1 (19) 
were found in the analysis of the charged lepton and quark mass matrices in refs. |pTl , |22 



The remaining ambiguity in the flavor independent quantity Iq can be removed by taking 
ratios. 

Note that the formulas in eqs. (|17D,(|18D are given in the field basis where the VEVs of 
the sneutrino fields are zero (2^1,2,3) = 0. Translation to this specific basis is achieved by the 



unitary rotation in the field subspace La = {Li, Hi) |1T7| , |30 |. 

The trilinear lepton couplings in eq. (|17|) must be antisymmetric under interchange 
of the first two indices i and j. This leads to the expression 

A.,.~^(e'-'°Af,-e'>'"»Ai). (20) 



The ratio of the Xijk and A^^^ is then given as 

Aijfc ( ^ _ h~u ^fk 

^ijk \^jk ^jk 



(21) 
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Approximating the ratio A^/ A:^ 



■yj^ — —2 (see ref [0) we get for the U{l)x charges given in 
eq. (p!9| ) the following couplings 



A' = a; 



333 



(22) 



A = Ao 



f (e^ - l)e^ \ 

(6-1)6 





, — A333. 



(23) 



V " / 

Here we denoted A' = {A^33},A = {Aj33}. In these eqs. remains only one free parameter 
A333. Thus, the U{l)x symmetry allows us to dramatically reduce the number of free pa- 
rameters in the neutrino mass matrix given by eqs. (^T])-(0). Totally, in the sector 
there are only four free parameters: the trilinear coupling A333 and the three bilinear A„ 
parameters (see eq. (pT])). The latter three are also subject to U{l)x constraints and under 
certain additional assumption further reduction of free parameters is possible |l7,30|. In our 
subsequent analysis we disregard these constraints and keep the three Ai^2,3 quantities as 
free parameters, avoiding additional assumptions. We already pointed out, that from the 
neutrino oscillation data alone one is able to fix the entries of the neutrino mass matrix up 
to the overall mass scale and the sign ambiguities which appear in solving the non-linear 
equations (|1),(|). As we have shown earlier in this section the MSSMxf/(l)x neutrino 
mass matrix depends only on four effective parameters. Now confronting the phenomeno- 
logical neutrino mass matrix Ai^'^, derived from the analysis of the neutrino data, with the 
MSSMxf/(l)x mass matrix A^'^, given by eqs. (p!T])-(p!^), we get a system of six linear 
independent equations 



(24) 



with five unknown quantities. Solving this system of equations one can uniquely determine 
the mass scale on the phenomenological side and in addition the absolute values of the four 
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theoretical parameters. The results for the neutrino masses with the corresponding CP- 
phases (cp and the family averaged Majorana neutrino mass (m^) = YliXi^ei''^i given 
in table 0. The ^ parameters are shown in table [m[ We made our numerical analysis 
under the assumptions discussed at the end of section |T|. Our results are given for the four 
different sets of matrix elements of the input matrix A4p^ found from the phenomenological 
analysis of the neutrino data in refs. [|lT|- [T^. Note, that for all the examined input sets 
we found a hierarchical neutrino mass scenario. 



The values of the bilinear parameters Ai 2,3 in table |T| are within the upper limits found 



in refs. |[T8|j3l|j3^ . The same is true for the values of the trilinear coupling A333 in table 
IT| as well as for all other couplings derived according to eqs. ( p2D and (pS)). They are not 
in conflict with the limits we found in section |T| and presented in table |. 

As seen from table H our prediction for the family averaged neutrino Majorana mass 
ranges in the interval | {m,y) \ ~ 0.01— O.OSeV. This range is more than one order of magnitude 



below the existing limit of |(?77.j,)| < 0.62 eV p^ , |25[] extracted from the current Oz//5/3-decay 
data [E3 . 



V. SUMMARY 

We discussed the phenomenology of neutrino oscillations and found upper limits on the 
entries of the three family neutrino mass matrix. This we did by using the double beta decay 
constraints on the family average Majorana neutrino mass (mjy). The so derived "maximal" 
neutrino mass matrix was used to test the flavor structure of the R-parity violating sector 
of the MSSM (R-parity violating Supersymmetric Standard Model). Comparing the 
theoretical 1-loop MSSM neutrino mass matrix with our phenomenological "maximal" 
matrix we extracted new limits on the trilinear ^p coupling constants Aj33 and A.33. These 
limits are more stringent than those existing in the literature. For the A233, A333 couplings 
are the new limits an improvement of up to 3 orders in magnitude compared to the existing 
limits. 
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As a next step we considered the ^pMSSM with the family dependent horizontal U{l)x 
symmetry. The framework of this MSSMx[/(l)x model is rather restrictive and allows 
one not only to set limits on the parameters from experimental data but, what is more 
interesting, to derive predictions testable in future experiments. Accepting the U{l)x charge 
assignment previously obtained in refs. PT| from the fit to the quark and charged lepton 



masses and mixing angles we related various trilinear coupling constants. As a result 
only four parameters in the MSSMxf/(l)x neutrino mass matrix remained free. In this 
case we were able not only to determine the upper limits but in addition the intervals of 
values for these parameters from the existing neutrino oscillation without any additional 
experimental information. Moreover, we completely reconstructed in this framework the 
three family neutrino mass matrix and give the predictions for the neutrino masses as well 
as for the family average neutrino Majorana mass (m,^). 

Noteworthy, our prediction for the (m^) lies in the interval ~ 0.01 — 0.05eV. This is one 
order of magnitude below the current experimental upper bound on this quantity. However 
the next generation 0^(3(3 experiments, e.g. the GENIUS experiment PB|, claim to be able 



to explore this region of small values of the average Majorana neutrino masses. Thus, we 
predict on the basis of the current neutrino oscillation data a positive result of searching 
for Oz//?/3-decay in the next generation experiments with the sensitivity to {m,y) in the range 
of few tens milli-eV. Since this prediction relies on the ^p MSSMxf/(l)x observation of 
OujSP-decay in this region would also indirectly support this model. 
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FIG. 1. The quark-squark (a) and the lepton-slepton (b) 1-loop contribution to the neutrino 
Majorana masses. The crosses on the hues denote the left-right mixing. 

TABLE I. Upper hmits for the trihnear R-parity violating couplings derived from the neutrino 
oscillation and the neutrinoless double beta decay data. Approximations are specified in section 

new limit Existing bounds (see 

MsusvllOOGeV ^-"^ ^ 3 -10 3 

Msus^fwOGeV '^^ ^ 6 -10 2 



A 



A; 



233 



3.8 •10"'^ 7 •10-'^ 

4.3 •10"'^ .36 



MsuSY/WOGeV 

MsusrflOOGeV ^.3 -IQ-^ .48 



TABLE IL The predictions of -^pMSSM with the U{l)x family symmetry for the neutrino masses 
rrii and the family average neutrino Majorana mass (my) = C^pC/g^mj. Here Cc^p,i = 1,2 are the 



CP phases of the neutrino mass eigenstates. The assumptions are specified in section [V. Different 



predictions correspond to different input sets for the neutrino mixing angles 6ij and Amfj found 
from the neutrino oscillation data in the papers cited in the last column. 

+ + .041 

+ .045 

+ + .010 HI 

+ .028 in 

+ .009 H 



mi[eV] 


7712 [eF] 


msleV] 


.004 


.032 


.549 


.018 


.036 


.549 


.002 


.002 


.030 


.000 


.0224 


.633 


.019 


.026 


1.054 
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TABLE III. The same as in table || but for the trilinear I^p coupling A333 and the bilinear 



parameters Aj defined in eq. (11) 



\^i\[GeV^ 
.008 
.008 
.004 
.006 
.004 



IA2|[GeF2 
.012 
.014 
.004 
.013 
.022 



|A3|[Gel/2 
.019 
.016 
.002 
.021 
.022 



-4| 



2.1 
2.4 
.7 
2.1 
3.0 



m 
m 
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